DNA Transport by a Micromachined Brownian Ratchet Device 
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We have micromachined a silicon-chip device that trans- 
ports DNA with a Brownian ratchet that rectifies the Brown- 
ian motion of microscopic particles. Transport properties for a 
DNA 50mer agree with theoretical predictions, and the DNA 
diffusion constant agrees with previous experiments. This 
type of micromachine could provide a generic pump or sepa- 
ration component for DNA or other charged species as part 
of a microscale lab-on-a-chip. A device with reduced feature 
size could produce a size-based separation of DNA molecules, 
with applications including the detection of single nucleotide 
polymorphisms. 



I. INTRODUCTION 

The Human Genome Project aims to provide the com- 
plete sequence of the 3 bilhon base-pairs of the human 
genome. While the dominant method for analyzing DNA 
fragments remains gel electrophoresis, new technologies 
that have the potential to increase the rate and decrease 
the cost of DNA sequencing and analysis, such as mass 
spectrometry and hybridization arrays, are critical to the 
project's success Q]. 

Here we describe a novel method of DNA transport and 
separation based on a Brownian ratchet. As described 
originally by Smoluchowski and noted by Feynman, a 
Brownian particle can undergo net transport on a poten- 
tial energy surface that is externally driven to fluctuate 
between several distinct states [||,^. Brownian ratchets 
have attracted theoretical attentio n [^ p^ due to their 
description of molecular motors p^l5[ and to their sim- 
ilarity with phenomena termed stochastic resonance and 
resonance activation |l6|Jl7tl . 

Brownian ratchets have been demonstrated to trans- 
port /xm- to mm-sized particles using dielectrophoresis 
psf , optical tweezers , and electrocapillary forces |2^] 
to generate ratchets. Other devices based on entropic 
ratchets or physical barriers |^,^ have been pro- 
posed as well. More recently, a geometrical sieve device 
has been used to separate phospholipids p^ . 

Despite these successes, the Brownian ratchet mecha- 
nism has not before proved capable of transporting DNA 
fragments in the size ranges applicable to DNA sequenc- 
ing (< 1000 nt) because the interactions used to estab- 
lish the ratchet potential were too weak. In contrast to 
previous devices using polarization interactions to gen- 
erate ratchets |l^Jl^, we have fabricated a device that 
uses charge-charge interactions to generate the ratchet 
potential. As seen below, the charge-charge interactions 
have sufficient strength to establish ratchets that can trap 
small DNA fragments. 

The ratchet-like wells that trap DNA are generated by 
charging a series of patterned electrodes. When the elec- 
trodes are discharged, the traps vanish and the molecules 
undergo Brownian motion. Next the potential is re- 
applied, and the particles again collect in the traps. A 
spatial asymmetry in the shape of each ratchet-shaped 
well rectifies the Brownian motion and produces net 
transport as the on-off cycle is repeated. Each molecule's 
transport rate depends on its diffusion constant, allowing 
the possibility of size-based separations. We have imple- 
mented the device by microfabrication on a silicon chip. 

This report describes the Brownian ratchet theory, 
provides a derivation of the transport rate, and presents 
experimental results for a single sized DNA oligomer. 
Greater details regarding the fabrication methods and a 



more extensive presentation of the experimental results 
for a variety devices and oligomer sizes are available else- 
where fE5|. 



II. THEORY AND METHODS 

A silicon wafer with six micromachined devices is 
shown in Fig. ^ , with a schematic design below. The 
electric potential that creates the ratcheting traps is gen- 
erated by two arrays of interdigitated electrodes that are 
perpendicular to the transport axis. The two sets of ar- 
rays each extend from bond pads on opposite sides of the 
device. The spacing between electrodes extending from 
the same bond pad is I. The asymmetric pattern creates 
two different spacings, r and I — r, between electrodes 
extending from opposite bond pads. The smaller spacing 
r defines the feature size of the device. 

As shown in also shown in Fig. ^, a simplified one- 
dimensional description of the potential approximates 
the electrodes as infinitely thin wires. To simplify the 
theoretical model, we have neglected the finite width of 
the electrode, the dependence of the electric potential on 
the distance normal to the surface, and the corresponding 
detailed calculation of the potential along the transport 
axis. 

When a voltage difference V is applied across the two 
electrode arrays during the on-state, with duration ion, 
sawtooth-shaped ratcheting traps are created for charged 
particles. The electrodes are discharged to V — during 
the off-state, with duration ioff, and particles undergo 
isotropic Brownian diffusion. When the potential is re- 
applied, the particles are again trapped in potential wells. 
The times ton and toff are within the low-frequency, quasi- 
static approximation; nonadiabatic effects and current 
reversals, reported elsewhere [^,^, are not applicable 
in this regime. Due to the asymmetry of the sawtooth 
shape and the choice of toe , a particle starting from well 
has a measurable probability to be trapped in well 1 
and virtually zero probability to be trapped in well — 1. 
Transport can be generated by repetitive cycling between 
the on-state and the off-state. The transport properties 
of a particle are determined by its diffusion constant D 
and its charge Q, along with the thermal energy fcsT and 
the device parameters defined previously. 

Each cycle begins by applying the potential V for a 
time ton that is sufficient to localize particles at the bot- 
toms of the trapping wells. A particle at the barrier-top 
drifts down the side of the sawtooth of length l — r to find 
the bottom of the well, which defines the time required 
for complete trapping. 
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according to overdamped Brownian motion. At the end 
of the trapping, the particle distribution at the bottom of 
a well is assumed to be much narrower than the feature 
size r. 

In the next phase of the cycle, the potential is turned 
off for time toff- When the potential is re-applied, parti- 
cles that have diffused further than the barrier to the left 
(roughly distance / away) will hop to the previous well, 
and particles that have diffused further than the barrier 
to the right (roughly distance r away) will be transported 
to the next well. Since r ^l, we can choose an intermedi- 
ate time toff such that toff /2D ^ P /2D and particles 
always move right, never left. 

To compute the probability a that a particle moves 
one well to the right in a single cycle, we note that the 
effective distance rcff it must travel is between r (the 
inner edge-to-edge distance between electrodes) and 3r 
(the outer edge-to-edge distance). The probability dis- 
tribution (or equivalently the Greens function) for a par- 
ticle starting at the origin undergoing one-dimensional 
Brownian motion for time toff is 



P(x;toff) = 



cxp[-xV4£'toff] 



off 



(2) 



An expression for a is then obtained by integrating 
P(a;;toff) from r^s to infinity: 



ierfc( 



rl^/4.Dtos) = ^erfc(v/tr/2toff). (3) 



The distance r^g has been written in terms of the char- 
acteristic diffusion time 2Dtr, where tr = r^ff/2D is the 
time required to diffuse a distance equal to the short 
side Toff of the trapping well. In this derivation, we 
assumed that toff is short enough that particles diffuse 
less than a single ratchet, i.e. P{x;tos) ~ for a; > L 
An expression for a valid in the limit of large toff is 

After each cycle, the particle distribution shifts to the 
right the distance al. After n cycles, the envelope of the 
distribution of particles in each well evolves as a Gaussian 
with center x{n) and square width cr^(n): 



x{n) = nla, 



tT^(n) = nPa{l — a). 



(4a) 



(4b) 



Here we have assumed that all the particles are in well 
at the start of the first cycle. 

Both a and the steady-state flux of particles through 
the device. 



flux 
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are plotted in Fig. ||. The transport fraction a (black 
line) increases with toff and approaches a maximum value 
of 1/2 when back-diffusion is neglected. The flux in units 
of the characteristic time tr is Q;tr/(ton+toff) and is shown 
for ton = ioff/3. The flux is non-monotonic and exhibits 
a maximum when toff ~ t^. Other voltage modulations 
more complicated than the on-off pattern described here 
are also possible and can change the direction of particle 
flow according to particle size [||. Other similar types 
of non-monotonic behavior have been termed stochastic 
resonance and resonance activation ||l6|,|l^ although they 
may also be described by dispersion and linear response 



Devices were fabricated from Pt, a non-reactive, 
corrosion-resistant metal chosen to avoid electrolysis of 
water, u sing relatively standard micromachining tech- 
nologies The fabrication began with thermally oxi- 
dized 100 mm diameter silicon wafers. A 200 A thick Ti 
layer was used as an adhesion layer between the subse- 
quent Pt layer and the silicon dioxide. Next, a 200 nm 
thick layer of Pt was deposited on top of the Ti. The 
electrodes were defined in the metal layers using pho- 
tolithography and ion milling. For the devices used in 
this work, the electrodes and the gaps between nearest 
electrodes were 2 /xm and the spatial period was 20/.tm. 

Oligomers labeled with fluorescent rhodamine dye 
(Amitof Biotech Inc., Boston, MA) were placed on the 
surface of the chip at 4 pmole//xl in deionized water. A 
microscope slide cover glass (Macalaster Bicknell, New 
Haven, CT) cut to size was used to confine the solution 
to a uniform thickness of approximately 10 /zm. A seal- 
ing compound was used to obtain a liquid tight seal along 
the edges, leaving the ends open. Square wave modula- 
tion (Synthesized Function Generator DS345, Stanford 
Research Systems, Sunnyvale, CA) with an amplitude of 
1.6 V and offset of 0.8 V (one set of electrodes at V 
and the second set at 1.6 V) was applied to the elec- 
trodes to generate the fiashing ratchet potential. Fre- 
quencies ranging from 0.7 Hz to 8 Hz were used, with a 
ratio ton/toS = 1/3. 

Video images of the analyte fluorescence were used to 
record the DNA transport. Images were captured us- 
ing a low fight imaging CCD camera (MTI VEIOOOSIT) 
mounted on an epi fluorescence microscope (Zeiss Ax- 
ioskop, Germany) using a lOx Fluor objective. The 
chip was illuminated using the output of a 50 W mer- 
cury lamp flltered with a green band pass filter. The 
brightfield fluorescence was imaged through a red low 
pass filter. The video images were recorded on video 
tape and transferred to a PC using a composite color 
PCI bus frame grabber (DT3153 Data Translation, Inc., 
Marlboro, MA). The fluorescence intensity resulting from 
the fluorescently labeled DNA fragments was analyzed in 



4 



a line across the video image synchronized to the ton pe- 
riod (HL Image-I— 1-97, Western Vision Software, Utah). 

According to Eq. |l], the expected time required to focus 
a DNA 50mer is 0.006 sec, based on a thermal energy of 
26 meV, a charge of l|e|~ per nucleotide yielding QV — 
80 eV, and a diffusion constant of 1.8 x lO^^cm^/sec (es- 
timated from a rhodamine-labeled 30mer on a quartz sur- 
face ) . We visually ascertained that the trapping time 
ton was sufficient to permit complete focusing of the DNA 
on the positive electrodes, even for the fastest switching 
rate of 8 Hz {ton = 0.03 sec). 



III. RESULTS AND DISCUSSION 

In Fig. ^ we show three images from a typical ex- 
periment using a device with 2 yum electrodes and a 
0.7 Hz switching frequency to transport a rhodamine- 
labeled DNA 50mer. These images were saved during the 
trapping phase of the cycle, and fluorescence maxima are 
clearly seen from DNA molecules captured on the posi- 
tive electrodes. At the start of the experiment, the DNA 
oligomers are focused on left-most three electrodes. As 
the potential cycles between on and off states, the packet 
moves to the right and broadens. The transport rate a 
can be estimated by noting that the intensity maximum 
moves from electrode 1 to electrode 3 after 10 cycles, then 
to electrode 5 after 10 more cycles, yielding a 0.2. For 
a more precise value, we extracted the intensity profile 
across the image, set a baseline at the 85**^ percentile of 
intensity, calculated the average position x{n) from the 
intensity above baseline, normalized x{n) by the 10-pixel 
spacing between electrodes, then measured the slope of 
x{n) to obtain a. For this experiment a — 0.18. 

The square width also increases linearly with the num- 
ber of cycles (data not shown); because we used a base- 
line threshold that narrows the width of the distribution, 
however, the formula of Eq. Eq no longer provides an 
accurate relationship between a^(n) and a. 

Fig. ^ shows the experimental results for a, along with 
±1(7 error bars from repeated runs. As predicted, a de- 
creases with increasing frequency. Also shown in Fig. ^ 
is a theoretical curve from Eq. g|. Least-squares fitting 
to the data from all but the highest frequency yielded 
D{r/rctf)^ — 3.5x 10~®cm^/sec. The good agreement be- 
tween the experimental results and the single-parameter 
theoretical fit supports our conclusion that transport is 
due to a Brownian ratchet. 

Furthermore, since r^s ~ (2-3)r, we find that D — 1.4- 
3 X lO^^cm^/sec. This indirect measure of the diffusion 
constant brackets the estimated diffusion constant for a 
DNA 50mer close to a glass surface, 1.8 x 10^''cm^/sec 

P9i. 



The experimental results demonstrating transport per- 
mit an examination of the feasibility fabricating a device 
to provide useful size-based separations of DNA. Separat- 
ing two chemical species requires that they have different 
diffusion constants D and D' and different hopping prob- 
abilities a and a'. The resolution between the species, 
defined as 



resolution = 



\x{n) — x'{n)\ 
0.5-[a{n) +a'in)y 



(6) 



improves as n^^^. The number of cycles required to reach 
the resolution of 1 typical for DNA separation applica- 
tions is 



a(l 



(a-a')2' 



(7) 



where we have assumed that the two packets acquire 
a similar width. The separation parameter tos (and, 
through toff, the quantities a and a') can be selected 
to optimize various quantities associated with a resolved 
separation, for example the device length, approximately 
I X [a/ {a — a')]^, or the total separation time, n x {ton + 

toff)- 

A typical application requiring the separation of DNA 
fragments is the analysis of single nucleotide polymor- 
phisms (SNPs). An SNP is a position in the genome 
where multiple nucleotides are likely. Characterizing ge- 
netic diversity through SNPs has applications including 
the identification of genes for disease inheritance and sus- 
ceptibility, the development of personalized medicines, 
and the documentation of human evolutionary history 
and migrations through a genetic record ||30|-|33|| . Prelim- 
inary sets of thousands of SNPs, identified by the 12 nt 
on either side of a polymorphism, have been reported 
p^ . Validating and detecting these SNPs can be ac- 
complished by resequencing specific 25 nt regions of the 
genome. 

Here we investigate the use of a Brownian ratchet de- 
vice for the resequencing and detection application. Cal- 
culations are based on a device with feature size r = 
0.1 /zm, periodicity 1 fim, and a potential difference of 
0.1 V between electrodes in the on-state. We estimate the 
effective DNA diffusion constant D{r/rcH)'^ from the the- 
oretical scaling for a self-avoiding walk, D ^ (length) 
p5| , ^ , and our experimental results for the 50mer. The 
persistence length of single-stranded DNA is 4 nm, or 
13.6 nt H^], indicating that the self-avoiding walk should 
be adequate (although not quantitative) for fragment 
sizes we consider. 

The detection of an SNP requires, at most, the abil- 
ity to sequence the 25 nt region surrounding the poly- 
morphism, which can be accomplished by separating a 
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DNA 24nier from a 25mer respectively. Using the the- 
oretical scaling of diffusion constant with DNA length, 
we extrapolate Dir/rcsY' of 5.44 x lO^^cm^/sec and 
5.31 X 10~^cm^/scc for the 24nier and 25iner. The time 
required to focus the DNA at the start of each cycle is 
ton ~ 1.8 X 10^'^sec. The optimized separation parame- 
ters, calculated using Eq. require toS = 2.5 x 10~''sec 
and 12,000 cycles for a total separation time of 5.4 sec 
on a 1.25 cm chip. 

In conclusion, we have fabricated a Brownian ratchet 
device that is capable of transporting small DNA 
molecules in aqueous solution, rather than the inconve- 
nient gel and polymer solutions required for electrophore- 
sis. This type of device could be used as a pump compo- 
nent for transport or separation of charged species in a 
microfabricated analysis chip. Multiple miniaturized de- 
vices can also be arrayed side-by-side for high-throughput 
operation. Based on experimental measurements, we 
suggest the feasibility of using this type of device for bio- 
logical applications, for example the validation of SNPs. 
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FIG. 2. The probability a that a particle hops one well to 
the right during a single cycle ol device operation is shown 
as a lunction of the ratio tr/tos (black line). For large toH, a 
approaches 0.5 as we do not consider back-diffusion. The flux 
in units of the characteristic time tr is atr/{ton + tott) and is 
shown for ton ~ tos/3 (grey line). The flux is reminiscent of 
a stochastic resonance with a maximum when tas ~ tr- 
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FIG. 1. The Brownian ratchet device is shown in 
schematic. Modulating the electric potential at the electrodes 
generates a ratchet-like potential energy surface for charged 
molecules like DNA. Cycling the ratchet between an on-state 
and an off-state generates transport. 
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FIG. 3. Three images are shown from a typical experiment 
using a device with 2 jj,m electrodes and a 0.7 Hz switch- 
ing frequency to transport a rhodamine-labeled DNA 50mer. 
These images were saved during the trapping phase of the 
cycle, and fluorescence maxima are clearly seen from DNA 
molecules captured on the positive electrodes. At the start of 
the experiment, the DNA oligomers are focused on left-most 
three electrodes. As the potential cycles between on and off 
states, the packet moves to the right and broadens. 
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FIG. 4. Experimental results for the transport of a DNA 
50mer by a device with 2 /im electrodes and 20 fim periodic- 
ity (points with la error bars) are compared with theoretical 
predictions (line). The theory requires a single adjustable 
parameter related to the diffusion constant of DNA. 
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